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1 Introduction 

Nonequilibrium kinetic Ising models , in which the steady state is produced 
by kinetic processes in connection with heat baths at different temperatures 
have been widely investigated [ 1 ] Most of these studies, however, are con- 
cerned with the effects the nonequilibrium nature of dynamics might exert 
on phase transitions driven by temperature. 

A different line of investigating nonequilibrium phase transitions has 
been via branching annihilating random walk (BAW) processes. Here par- 
ticles chosen at random carry out random walk (with probability p) with 
annihilation upon meeting. The increase of particles is ensured through 
production of n offsprings with probability 1 — p. It has been shown by 
Sudbury [2] , that the n = 2 BAW in one dimension leads to extinction for 
any p > 0. Otherwise a phase transition occurs for finite p, which is in the 
same universality class as directed percolation (DP) if n = odd, while the 
critical behavior is different for n = even[3 ]. 

Grassberger et al. [ 4 ] studied probabilistic cellular automata models in 
one dimension involving the processes k ^ 3k and 2k —^ {k stands for 
kink), very similar to BAW with n = 2. These models, however, do show 
a phase transition and both time-dependent and steady-state simulations 
have resulted in non-DP values for the relevant critical exponents [5 ]. 
Quite recently, Jensen [6 ] has reported computer simulation data according 
to which the n = 4 BAW is in the same dynamic universality class as 
Grassberger 's cellular automata. 

It is the aim of the present paper to introduce a class of general nonequi- 
librium kinetic Ising models (NEKIM) with combined spin flip dynamics 
at T = and spin exchange dynamics at T = oo in which, for a range 
of parameters of the model, Grassberger-type transition takes place. The 
advantage over the cellular automaton formulation is that in NEKIM the 
rates of random walk, annihilation and kink-production processes can sep- 
arately be controlled. This circumstance leads to the understanding of the 
n = 2 BAW result which occurs in our model as a special case. Results 
of computer simulations are presented for different critical exponents and 
scaling function from random as well as 1-kink initial states. 



2 The model 

The model we will investigate here is a one-dimensional kinetic Ising model 
evolving by a combined spin-flip and spin-exchange dynamics with the re- 
striction that spin flip transition probabilities satisfy detailed balance for 
the equilibrium state of the (ferromagnetic) Ising model at temperature 
T = while the spin exchanges are random (T = oo). 

A general form of the Glauber spin-flip transition rate in one-dimension 
for spin Sj sitting at site i is [7] (sj = ±1): 

Wi = —(1 + (5Si_lSi+l) (l - -Si{Si_i + Si+i) j (1) 

where 7 = tanh2J/A;T , with J denoting the coupling constant in the Ising 
Hamiltonian F and 6 are further parameters which can also depend on 
temperature. There are three independent rates: 

Wsan^e = ^(l + 5)(l-7) (2) 

^oppo = ^(l + 5)(l + 7) (3) 

Windif = ^(^~^)' (^) 

where the suffices same etc. indicate the three possible neighborhoods of 
a given spin(t|t, ||J, and t||, respectively). In the following T = will be 
taken, thus 7 = 1, Wsame = and F, 6 are constants to be varied. In this 
limit the three best-known single spin flip kinetic Ising models correspond 
to the following choices of parameters : 
a). Glauber model [8 ]: F = 1, 5 = 
b). Metropolis model[9 ]: F = 2/3, 6 = -1/3 
c). Haake-Thol model[10 ]: F = 2, 5 = 1 . 

The spin-exchange transition rate of neighbouring spins, originally in- 
troduced by Kawasaki [7 ] to situations with conserved magnetization, can 
be given, for 7 = 0, as: 

Wii+l = -Pex[l - SiSi+i] (5) 

where Pex is the probability of spin exchange. 



The transition probabilities (2)-(4) are responsible for basic elementary 
processes of kinks in the usually expected ferromagnetical ordering of Ising 
spins. A kink separating two domains can carry out random walks with 
probability 

Pr^ oc 2wi„di/ = -(1-5) (6) 

while two kinks getting into neighbouring positions will annihilate with 
probability 

Pan OC Woppo = -{I + S) (7) 

{wgame IS responsible for creation of kink pairs inside of ordered domains at 

In case of the spin exchanges, which also act only at phase boundaries, 
the process of main importance here is that a kink can produce two off- 
springs by the next time step with probability 

Pfc-.3fc OC Pex (8) 

The abovementioned three processes compete, and it depends on the pa- 
rameters r, 6 and Pex what the result of this competition will be. It is 
important to realize that the process k ^> ?>k can develop into propagating 
offspring production only Mprw > Pan, i-e. the new kinks are able to travel 
on the average some lattice points away from their place of birth an can 
thus avoid immediate annihilation. It is seen from the above definitions 
that (5 < is necessary for this to happen. In the opposite case the only 
effect of the k ^ 3k process on the usual Ising kinetics is to soften domain 
walls. As a matter of fact, in earlier investigations (though from a different 
point of view) of the competing Ising kinetics in one dimension the choices 
of parameters correspond to the latter situation [11]. 

3 Phase boundary, critical exponents 

We have considered a simplified version of the above model by keeping only 
two parameters instead of three by imposing the condition 

Pk^Sk = 1 - {Prw + Pan)- (9) 



In the plane of parameters Pex and 1 — 6 = -.i''''" ^ " the phase diagram 

^~rPrw /Pan 

shown on Fig.l. has been obtained by computer simulation. The initial 
state has been random with zero average magnetization . Application of 
our nonequilibrium rule then corresponds to quenching from T = oo to 
T = with subsequent exchanges (T = oo) at each time step. The line 
of phase transitions separates two kinds of steady states reached by the 
system for large times: in the Ising phase the system orders while the 
active phase is disordered from the point of view of the underlying spins. 
The cause of disorder is the steadily growing number of kinks with time. 
Kink-multiplication is here a branching process, as explained above, in a 
sense similar to directed percolation while in the usual case of Ising phase 
transition with Glauber kinetics kinks are created pairwise inside of ordered 
domains at nonzero temperature. 

Exponent a, used for determining the phase boundary, as well as other 
exponents occuring below are defined in the framework of the scaling con- 
siderations presented by Grassberger [5] . The branching increase of kinks 
becomes obvious when starting from a single kink state at t = 0. The 
scaling form for the density p(x, t, e) of the kinks has been supposed in the 
form [5] 

p{x, t, e) oc t-"0(ex^/'^^ , et^/^n ) (10) 

Here e measures the deviation from the critical probability at which the 
branching transition occurs, u^ and z/y are exponents of coherence lengths 
in space and time directions, respectively. 0(a, b) is analytic near a = and 
6 = 0. Using (10) the following relations can be deduced. 

The average particle number N{t, e) grows in the active phase as 

A^(t,e)oct''^(eti/"ii). (11) 

The RMS size of the cluster growing from the single-kink initial config- 
uration is given by 

1 /9 

< (X- < X{t) >f> oc < Xmaxit) - Xmin{t) > OC t^^^ (12) 

where Z = 2z/_|_/z/||. 

When starting from a random initial state the exponent f3 characterizes 
the growth of the average kink density in the active phase: 

p(e) = Um p{x, t, e) oc e^ (13) 

t — *oo 



while the decrease of density at the critical point is given by 

n(t) oc t~". (14) 

The exponents are connected by scaling laws: P = u^a, a = Z/2 — rj. 

The phase boundary shown on Figl. has been identified using (14). It 
is worth noting that the critical point marked on the phase boundary has 
been chosen in a region where the width of the critical region is relatively 
small. The rest of the critical exponents have been measured only around 
this point :^^ = .35, 4 = -.4 ± .01 and (pex)^ = 1 - ft = -3 ± -Ol- The 
deviation from the critical point has been chosen in the direction of 6. 

The critical region is especially wide near the two ends of the curve. 
Moreover, near pex = 1-0 where we get close to 5 = (Glauber case), it 
is very hard to determine a: the exponent grows slowly with time from 
a nearly zero value at small times. Several runs with different values of 
the parameters also without the restriction (9), have been performed in 
this Glauber limit all showing that the 6 = case for all F, pex values, 
remains Ising-like: the exponent a tends to the value a = .5 for large 
times. The accuracy of our simulations does not allow us to make a quite 
definite statement, nevertheless it seems very likely that the asymptote of 
the phase boundary is 5 = for pex = 1-0. Therefore the conjecture that 
Prw > Pan IS a uecessary condition for a BAW-type phase transition to occur 
seems to be supported by the present simulations. 

The value of a provided by our simulations agrees, within error, with 
that obtained by Grassberger[4,5] in a one -dimensional elementary cellular 
automaton model with special added noise p. The same is true also for the 
rest of the critical exponents reported below. Fig. 2 shows the density of 
kinks in the active phase for different values of the deviation from the 
critical point e =| S — 6c |, starting from random initial states. At e = the 
power law behaviour (14), is seen with a = .27 ± .02. [3 has been obtained 
directly, according to eq.(13), by taking the level-off values of n{t), some 
of which can be read off from Fig. 2. Though this is known to be a fairly 
inaccurate way to determine (3 [5], we have found -though only over one 
decade in e - the reasonable value P = .8 ± .08. 

On Figs. 3 through 5 results of simulations starting from one-kink initial 
states are seen. Data for N(t, e) at e = are presented on Fig. 3; the straight 
line on the log- log plot leads to ?] = .30 ± .02. Fig. 4 shows the growth of the 
average cluster size, eq.(lO), again on log-log scale, yielding Z/2 = .57±.01. 
Supposing the scaling law j3 = u\\a to be valid we get u\\ = 2.9. 
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Using the values obtained for the exponents, we have checked the con- 
sistency of the data within the scahng framework by measuring the scahng 
form, eq.(ll). The result is seen on Fig. 5 for some values of e, nevertheless 
also only inside of a decade, because of computational difficulties. The sec- 
ond scaling law -.a = Z/2 — rj is also fulfilled with the exponents found in 
our simulations. 



4 Connection with BAW 

The definition of BAW is the following. The process is taking place in 
two steps: a), a particle chosen at random moves to a nearest neighbour 
place also randomly chosen, with probability p and b). a particle chosen 
at random produces n offsprings on neighbouring sites with probability 
1 — p. If any of the neighbouring places in question is already occupied 
then annihilation takes place. The k + k ^ process is prescribed to 
occur with probability unity. As mentioned in the introduction, no phase 
transition takes place in the n = 2 case. 

We can bring into relation the n = 2 BAW in Id and NEKIM by noticing 
that in the BAW rule the transition probability for annihilation is unity. 
From eq.(7) Pan = 1 leads to F = (1 + 5) giving, with eq.(6), Prw = jr^- 
From the trivial correspondence 

Prw = P (15) 

Pk^sk = l-p (16) 

while changing p, the parameters of the corresponding NEKIM models take 
the values from 6 = 1, pex = 1, F = 2atp = through values 5 > 0, pex < 1, 
F < 2 to (5 = 0, Pex = 0,F = latp = l. In these regions, however, no phase 
transition has been found to take place, all steady states are Ising-like. 



5 Discussion 

Non-thermal phase transition has been found in a family of one- dimensional 
kinetic Ising models evolving by a combined spin-fiip and spin-exchange 
dynamics. The single spin fiip transition rate is supposed in its most general 
form. It contains three parameters of which the temperature is fixed at 
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zero where the remaining two control the probabiUty of random walk and 
annihilation. The nonzero transition rate of spin exchange leads to A; — > 3A; - 
type kink production and provided RW overwhelms annihilation, increasing 
rate of kink production results in a phase transition from Ising-like to active 
phase. Numerical simulations have given evidence that the model is in 
the same universality class as Grassberger's cellular automata models[4,5], 
where a certain stochastic element (with probability p) of a simple Id CA 
rule drives the transition. For p = the rule is deterministic and the steady 
state is infinitely degenerate, when starting from a random initial situation 
the system settles in one of these states after a few iteration steps. The 
corresponding kinetic Ising model is the Haake-Thol model in the T — > 
limit mentioned in section 2. Increasing p up to Pc in the CA model, could 
correspond to a line on our phase diagram (Fig 1) starting at the origin 
and ending at some intermediate point of the transition line. It would be 
difficult, however, to make a closer (quantitative) correspondence between 
the two kinds of models since in CA this single parameter p is responsible 
for random walk as well as A; — * 3A; in a complicated, nonlinear way. To find 
correspondence between our model and the two offspring BAW is more 
straightforward. For this aim restriction (9) had to be lifted since pan = 1 
is taken in BAW rules and then eq (9) would lead to negative probabilities. 
According to the correspondence between BAW and NEKIM it is clear that 
a necessary condition to find a phase transition is Pan < 1- 
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Figure captions 

Fig 1 Phase diagram of the two-parameter model. The phase boundary 
has been obtained by measuring n{t), the density of kinks, starting from 
a random initial distribution and locating the phase transition points by 
n{t) oc t~° with a = .27 ± .04. Typically the number of lattice points 
has been / = 2000 and averaging over 500 independent runs has been per- 
formed. Checks for unimportance of finite size effects have been carried 
out. The dotted line with arrow indicates the critical point and the direc- 
tion of deviation from it along which further critical characteristics have 
been determined. 

Fig 2 Density of kinks as a function of time on a log-log plot for dif- 
ferent values of e =| S — 6c |, where 6c = —.4 ± .01 is the critical value 
of parameter 6 at the chosen critical point. Pex is kept fixed at its criti- 
cal value. Lattice sizes : 2000—8000, number of independent runs :400— 800. 

Fig 3 The growth of the average particle number at the critical point 
when starting from a single-kink initial state. The number of independent 
runs in the averaging: 900. The scale is double logarithmic. 

Fig 4 Average distance between rightmost and leftmost kinks at the criti- 
cal point on a double logarithmic scale. Average over 900 independent runs 
has been carried out. 

Fig 5 Scaling function ^(et^'^'n) of eq.(ll) for different values of e. Ini- 
tial state: single-kink, number of independent runs in averaging : 500. 
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